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Abstract 

We consider the problem of hypotheses testing with the basic sim- 
ple hypothesis: observed sequence of points corresponds to stationary 
Poisson process with known intensity. The alternatives are stationary 
self-exciting point processes. We consider one-sided parametric and 
one-sided nonparametric composite alternatives and construct locally 
asymptotically uniformly most powerful tests. The results of numeri- 
cal simulations of the tests are presented. 
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1 Introduction 

Let {^1,^2, . . .} be a sequence of events of a stationary point process X = 
{Xt, t > 0} {Xt is a counting process). The simplest stationary point process 
is, of course, Poisson process with a constant intensity 5* > 0, i.e., the incre- 
ments of X on disjoint intervals are independent and distributed according 
to Poisson law 

Qk (f _ ^]'' 

P{Xt-X, = k}= \ e-^^'~'\ 0<s<t, A; = 0,1,... 

Therefore if we have a stationary sequence of events it is interesting to check 
first of all if this model (Poisson process) corresponds well to the observations. 
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The importance of this problem was discussed by Cox and Lewis (1966), 
Section 6.3. 

The alternatives close to the basic hypothesis correspond to the case when 
the non-poissonian behavior is due to small perturbations of the Poisson 
process and are the most interesting to test. For "far alternatives" any 
reasonable test has power function close to 1 and the comparison of tests 
seems less important. Let us consider the problem of small signals detection 
by the tests of fixed size e e (0,1). Using the terminology of statistical 
radiothechnics we say that there is at least two types of close alternatives: 
the first one corresponds to small "signal-noise ratio" (signals of small energy) 
and the second, when the amplitude of the signal can be small, but the total 
energy due to the sufficiently long time of observation is comparable with the 
noise energy (see, e.g., Kutoyants (1976)). For the first class of alternatives 
the approach of locally optimal tests, which provides the optimality of the 
power function at the small vicinity of the basic hypothesis (the values of 
the power function are close to e) was developed (see, e.g. Capon (1961)) 
and for the second class of contiguous alternatives the optimality of the test 
for a wider class of close alternatives (the values of the power function are in 
(s, 1)) is proved (Pitman's (1948) approach, Le Cam's (1956) theory). 

For stationary point processes with Poisson hypothesis and stationary 
alternatives Davies (1977) proposed the locally optimal (efEcient) or asymp- 
totically locally efRcient test. This test is based on the comparison of the 
derivative of the log-hkelihood ratio with some threshold. See as well Da- 
ley and Vere- Jones (2003), Section 13.1, where the approach of Davies is 
discussed. 

In the present note we suppose that we have observations of the point pro- 
cess X'^ — {Xt, < t < T} on the interval [0, T] and consider two problems 
of hypotheses testing in the asymptotics of large samples (T — > oo). In both 
problems the basic hypothesis is simple: the observed process is standard 
Poisson with known constant intensity 5* > 0. The composite alternatives 
are: the observed process is a realization of self-exciting point process (some- 
times called Hawkes (1972) process) with in the first case intensity function 
depending on one-dimensional parameter and in the second case the intensity 
function belonging to a wider (nonparametric) class of functions. We follow 
the mentioned above Pitman-Le Cam's approach. We start with the locally 
asymptotically uniformly most powerful test (LAUMPT) in the parametric 
case and the main result of the presented work is the LAUMPT where the 
optimality is shown for sufficiently large class of local nonparametric alter- 
natives. The similar results for diffusion processes can be found in lacus, 
Kutoyants (2001) (small noise asymptotics) and Kutoyants (2003) (ergodic 
processes). 
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2 Preliminaries 



Remind several facts from the theory of point processes (the details can 
be found, for example, in Liptser and Shiryaev (2001), Chapter 18). Let 
{fl, ^, P) be a probability space and let {^t, t > 0} be a nondecreasing family 
of right continuous a-algebras G d for any < s < t. We denote by 
ti,t2, . . . a sequence of Markov stopping times adapted to {^t, ^ > 0} (that 
means {u : ti < t} & for all t > 0). Let Xt be the number of events ti up 
to time t, i.e., X = {Xt,^t,t > 0} is a random process such that 



where y, is the indicator-function of the event A. 

^{A} 

We assume that EX^ < oo (there is no accumulation points on any 
bounded interval). The process X admits a unique (up to stochastic equiv- 
alence) decomposition (Doob-Meyer decomposition) 

Xt = At + Mt, (1) 

where Ai = {M-tydtyt > 0} is a martingale and A = {At,'^t,t > 0} is pre- 
dictable increasing process (Liptser and Shiryaeyv (2001), Theorem 18.1). 
We suppose that the compensator A is absolutely continuous 

At= [ S{v,uj) dv, t>0 
Jo 

where S = {S {t,uj) > 0} is called intensity function. We suppose as 
well that ([T]) is the minimal representation of the point process, i.e., S (t, uj) 
is measurable w.r.t. a-algebra generated by {Xs, s < t} for any t > and 
we write S(t,u) = S{t,X). To describe a point process it is sufficient to 
specify its intensity function. We study in this work a special class of point 
processes with intensity functions which can be written as stochastic integrals 
with respect to the past of the underlying point process. 

In the particular case when S is deterministic, the process X is (inho- 
mogeneous) Poisson process with intensity function S{v,uj) = S{v). In this 
case 



/; s (v) dv 



k 



I>{Xt-X, = k} = ^ ^exp<^-/ S{v) dv 



k\ 

for any t > s > and A; = 0, 1, . . .. If the assumption of the independence of 
increments is no more valid, then S is no more deterministic and X can be 
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a stationary point process (see Brillinger (1975) and Daley and Vere- Jones 
(2003) and references therein for wide classes of such processes and their 
applications in real problems). 

(T) 

Remind that the distribution of the point process in the space of its 
realizations (^(0,T) , ^t) is entirely characterized by its intensity function 
S. The likelihood ratio formula (w.r.t. Poisson process of constant intensity 
S^) has the following form (see Liptser and Shiriyev (2001), Theorem 19.10) 

L (X^) = exp In dX^ - £ [S (i, to) - S.] di | , 

where we suppose that the intensity S (tjCu) is left continuous function and 

S {t,uj)dt <oo^ 

under all alternatives studied in this work. 



3 One-sided parametric alternative 

Suppose that we observe a trajectory X'^ — {Xf, <t <T} oi point pro- 
cess of intensity function St (i?) ^ {S {'&,t,uj) ,0 < t < T}. If = 0, then 
S (0,t,cj) = 5*, i.e., this point process is Poisson process of intensity S** > 0. 
Under alternative > and St {'&) is the intensity function of self-exciting 
point process. As usual in such problems, we consider contiguous alternatives 
(Pitman's (1948) alternatives, Roussas (1972)), hence we change the variable 
-& — u/VT and test the following two hypotheses 

Jifo: u = 

jTi : -u > 0. 

We denote Eq the mathematical expectation under the hypothesis J^, 
and E„ under (simple) alternative -d = u/\^. 

Let us fix £ G (0, 1) and denote by the class of test functions (pT {X^) 
of asymptotic size e, i.e., for (pT £ we have 

lim Eo (pT (X^) = e. (2) 

As usual, (pT {X^) is the probabihty to accept the hypothesis having 
observations X^. The corresponding power function is 

(3j, {u, (Pj) = E„ (Pt {X^) , u>0. 
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We introduce the asymptotic optimality of tests with the help of the 
following definition Le Cam (1956). 

Definition 1. A test 0^ (■) is called locally asymptotically uniformly most 
powerful in the class if for any other test (•) € and any constant 
K > we have 

lim inf T/?^ (u, 6*) — Brj. (u, > 0. 

Our goal is to construct locally asymptotically uniformly most powerful 
test in class Jt^. 

Self-exciting type processes were introduced by Hawkes (1972) and de- 
fined by intensity function of the following form 

S{t,u;)^S,+ f g{t-s) dX,^S, + J29{t-ti), (3) 

where 5"* > 0,ti are the events of the point process and the function g{-) > 
satisfies the condition 

poo 

p=/ g{t)dt<l. (4) 

Remind that according to this representation of the intensity function, the 
distribution of ti is exponential at rate and for all n > 1 

P {tn+i > t \ ti, . . . ,tn} ^ exp ^-SJ -J Y^9{s-ti) dsj . 

Note that A (t) — FiXt is solution of the equation 

rt pt pv— 

A{t) ^ E S {v,uj) dv ^ SJ + E g{v- s) dX, dv = 

Jo Jo Jo 

^ SJ+ f f g{v-s)A'{s) dvds. 
Jo Jo 

In stationary case the intensity S {t, a;) , is a stationary process 
S{t,u;) =S^+ f g{t-s) dX„ 

J —oo 

and 

A (t) = t = fjLt. 

1- p 
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The spectral density of this process is 

/(A) ^ 



27r|l-G(A)r 
where 

/•oo 

G{X)= / e'^'git) dt, p = G(0). 

Example 1. Let g{t) — ae"'''*, where a > 0,7 > and a/^ < 1. Then the 
point process X with intensity function 



S {t, u) = S^, + a ^ 
is self-exciting with the rate 



7 — 0; 



Example 2. The function g[-) can be chosen in such a way that the spectral 
density of the point process will be rational 

|g(iA)|^ 

where Q {z) = z'f + aiz^'^ + . . . + and P {z) = z'f + hiz^'^ + . . . + 6p. It is 
supposed that P {■) and Q{ ) have no zeroes in common and no zeroes in the 
closed right half plane (see Pham (1981), where the asymptotic properties of 
the MLE for this model are described) . 

We assume that the observed process is either Poisson with constant 
intensity or self-exciting with contiguous intensity function 

S{^,t,u})^ S^ + '&T [ h{t-s) dXs. 
Jo 

Contiguous means, that the likelihood ratio is asymptotically non degenerate. 
The function h{-) is supposed to be known, bounded and 

/i(-)e£^(M+) = |/(-)>0 : J^^ fit) dt<oo|. 
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To have contiguous alternatives we choose, as usual in regular problems, 
§T = u/ -\/T, i.e., 

S{u,t,uj) = S^^ — -1= / h{t — s) dXs, u > 0. 
vT Jo 

Note that for any h (■) G C\. and any u < K for sufficiently large T the 
condition (jlj) is fulfilled for the corresponding function g {■) = uT~^^'^h (■). 
This leads us to the following one sided hj^otheses testing problem: 

: u = 0, (Poisson process) 
J^i : M > 0, (self-exciting process) . 

This model corresponds to small self-exciting perturbations of the Poisson 
process of intensity S^:. 

Note that as we use the LAN approach, we study the behavior of the tests 
statistics under hypothesis only (Poisson process with constant intensity) and 
do not use the stationarity of the self-exciting processes under alternatives. 
The limit of the power function is obtained using LAN and Le Cam's Third 
Lemma. 

Let us denote 

At (X"^) = r f hit- s) dX, [dX, - dt] . 



Here 



[ h{t - s) dXs = ^h{t - ti) 
Jo . w 



ti<t 

(limit from the left of the integral, i.e., the term with Sj = t is excluded) and 

2 



h{tf dt + S, U h{t) dt 



is the Fisher information of the problem. Throughout this paper we denote 
by the 1 — e quantile of the Gaussian law A/'(0, 1). 

Theorem 1 Let h{-) E C\_ (M_|_) and bounded. Then the test 



K (^'^) - ^{A^(xr)>.,} 



is locally asymptotically uniformly most powerful in the class and for any 
u > its power function 



where C ~ A/'(0, 1). 
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Proof. First note that the family of measures |p^"^'','i9 > o| under hypoth- 
esis ^0 is LAN at the point = 0, i.e., the random function Zt{u) = 
L (u/y/T, X^j admits the representation (see Kutoyants (1984), Theorem 



4.5.3) 



Zt{u) = exp<j 1^ + h{t-s)dX,] dX,- 



exp |« v^A^ (X^) -^ll + r^ («, X^) | 



where 



Co{A^{X^)}^Af{0,l) (6) 



and Trp (^Uj,,X^) — >• for any bounded sequence {m^}. 
To verify we check the following two conditions: 

• Lindeberg condition for stochastic integral: for any S > 

1 



the law of large numbers: 



1 



Here we denoted 



Ht= h{t-s) dX,. 

^0 

By these conditions the stochastic integral At (X-^) is asymptotically nor- 
mal. The proof of the corresponding central limit theorem can be found, say, 
in Kutoyants (1984), Theorem 4.5.4 (of course, this theorem is a particular 
case of general CLT for martingales). 

To check these conditions we introduce an independent Poisson process 
{Xt, t < 0} of intensity and replace Ht by 

= f hit- s) dX,. 



It is easy to see that for the process Hl,t > we have 



Eoi/; = / h{v) Av 
'o 



and 



Eo ( [H; - Eo HI] [H: -EoH:]) = S, h{v + t-s) h (v) dv. 

J max(0,s— t) 

Note as well that 







<^ f Eo( [ h{t-s) dX^ dt 



Tz/2 





2 /-T 



oo / foo \ 2 

/i (f )^ dv + i h (f ) dt> 



dt — >0 



as T ^ oo. 

Now the process Hl,t > is second order stationary and 

/oo / POO \ 2 

/i (t)' dt + sUj h (t) dt j = Eo {Hlf < oo. 



Hence 



as T — > oo and 

JL™J^/E„(i/-X{,„,,,^,) d( = 0. 

The law of large numbers ([7]) will follow from the convergence: 

Mt = Eo(^^£ Hfdt-Eo{H*f^ = 

= ^2^ £ Eo {Hf - Eo {H*f) {Hf - Eo (//o*)') dt d. 
To prove it we need the following elementary result. 



Lemma 1 Let X = {Xf, t G A} he a Poisson process of constant intensity 
S>QonA^R,andletf{-),g{-)eC''{A) = ^f{-) : j \f {t)\^ dt < oo^ , 
A; = 1,...,4. Then 



Covl ( / f{v) dXA ,( [ g{v) dX„' 



= 4 /" f{v)Sdv I g{v)Sdv [ f{v)g{v)Sdv+ 

J A J A J A 

+ 2 (^^ f{v) g{v) Sdv^ f{v) g\v) S dv+ 

+ 2 / f{v)Sdv f f{v)g'{v)Sdv + 2 f g{v) S dv [ f{v)g{v)Sdv. 

J A J A J A J A 

Proof. Using well-known properties of the Poisson processes (see, e.g., Ku- 
toyants (1998), Lemma 1.1), we obtain the moment generating function 

(A, ii) = Eo exp j-^ y / {v) + J 9{v) dX^ 
= exp jy (e^/(^)+''ff(^) - 1) 5 dvj . 



Remind that 
Gov 



fiv) dXj ,( [ g{v) dX,^ 



(A, fi) 



(A, 0) 



A=0,/i=0 



(0, H) 



A=0 



M=0 



Therefore the proof of the lemma follows from direct calculations. 

Now we can write 

R (t, s) = Eo {{H:f - Eo {H*f) {{H:f - Eo {H^f) = 

/tAs 
h{t-vfh{s-vfdv+ 
-oo 



tAs 



[h {t - vf h{s - v) + h{t - v)h{s - vf] dv 
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where we put 

a = I h{y) dy 

and (for t = t — s) 







POO 

Kit,s) = S, / hiy) hiy-\T\) dy = K (r) . 

J\r\ 

Further, as the function h{-) is bounded, we have the estimate 

R{t,s) <C K{t). 

Hence 



Mt = ^ j\{t,s) dtds<^ £ j\{t,s) dtds< 



T 

For the function K (■) we have 



< - Kir) dr. 



-T 



/T i>T POO / POO \ 2 

K{T)dT = sJ / h{y) h{y-\T\) dydT<2sA h{y) dy] . 
-T J-tJ\t\ \Jo 



Hence Mt — > and we have the law of large numbers (j7]) . 

The property 0t (■) ^ follows from the mentioned above asymptotic 
normality of the statistic {-^'^) ■ 

Note as well that the convergence follows from 

(see the Third Lemma of Le Cam (van der Vaart (1998), p. 90)). 

The asymptotic optimality of the test follows as well from the general 
theory (see, e.g., Le Cam (1956) or Roussas (1972)), because if we replace 
J^i by any simple alternative ^ : u = u^:, then the test 

is the most powerful. Here 

b, = exp - ^ul I* I (1 + (1)) . 
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It is easy to see that (pT (■) £ ^-nd the power function 

Therefore the test 0t (•) is asymptotically as good as the likehhood ratio test 
for any simple alternative. 

Remark 1. Note that the statistic (-^^) can be written as follows 



ho<ti<Tti<tj V h Q<t^<T''^ 



h {v) dv, 



where U are the events of the observed process. 



Remark 2. By a similar way we can consider the problem of contiguous 
hypotheses testing when under the hypothesis Mq the observed process is 
self-exciting too. For example, let h{'d,t) > 0, i > be a smooth function of 
such that for all ■»? e © the condition 



/ h {'d, t) dt<l 
Jo 



holds. Then with the help of this function we introduce a family of self- 
exciting processes with intensity functions 

S{'d,t,u})^S^+ [ h{'d,t-s) dXs. 



Remind that these are stationary processes. 
Now we can test the hypotheses 

by the observations = {Xt, <t <T}. Suppose as well that the function 
h •) is two times continuously differentiable on i) at the point i) = 'do and 
the derivatives /i (i?, •) , h{'t}, ■) satisfy the suitable conditions of integrability. 
Let us denote 

(^) = £ h{^,t- s) dX,, (^) = £ dX,, 

and put 

^ £ srnk) - - ^' ^'°> ■ 
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Then it can be easily shown that the test 

where = Z£^/Th{^^ is chosen from the condition 0r G is locaUy 
asymptoticaUy uniformly most powerful in the class Here Ih (i^o) is the 
Fisher information 

where ^ {'do) and ^ (^9o) are stationary random variables related to the limit 
distribution of the vector (''^o) , Ct (''^o)- 



4 Testing of dependence 

Suppose that we have two sequences of events < ti < t2 < ■ ■ ■ < < T and 
< Si < S2 < ■ ■ ■ < sm < ^ with corresponding counting processes X^ — 
{Xt, < t < T} and = {Yt, 0<t<T}. The first process is Poisson with 
constant known intensity function Sx {t, uj) = Sx > and the intensity 
function of the second process can be written as 

Sy {t, uj) = Sy+ j r{t-s) dX„ 

where r (•) G (1^+)- Therefore, if r{t) = 0, then the observed processes 
are standard (independent) Poisson processes of intensities Sx and Sy re- 
spectively (Hypothesis Mq). For the other values of r (•) we have dependent 
point processes. 

We suppose that the dependence between these two processes, if exists, 
is weak, i.e., the function r (•) is sufficiently small and we can apply the local 
approach. As before we suppose that r {t) — -dr h {t), where h (•) e {^+) 
and '&T = u/VT — >• 0. 

J^fo '■ u — 0, (independent Poisson processes) 
Ml : u> 0, (depending processes) . 

Introduce the statistic 

{X^, Y^) = f f h{t- s) dX, [dY, - Sy dt] = 

Jo Jo 
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where 

Q / poo / /'OO \ 2\ 

Proposition 1 Let h (•) e C\, and bounded. Then the test 

is locally asymptotically uniformly most powerful in the class and for any 
u > its power function 

^/3(m) = p|C>2;.-m^}, (8) 

where C ~ A^(0, 1). 

Proof. The proof is quite close to the given above proof of the Theorem 1, 
and hence is omitted. 

Remark 3. The similar problem can be considered for the couple of mutually 
exciting point processes with intensity functions 

Sx {t, u)^Sx+ [ rxY {t - s) dY„ 

J —oo 

Sy (t, u;)^Sy+ f rYx (t - s) dX„ 

J —oo 

where rxY (■) j^yx (■) £ Therefore, if Vxy (t) = and ryx (t) = 0, 

then the observed processes are standard (independent) Poisson processes 
of intensities Sx > and Sy > respectively (Hypothesis Jifo). Under 
alternative there exists a weak dependence of these processes through their 
intensity functions. 

5 One-sided nonparametric alternative 

In all considered above problems the alternatives are one-sided parametric. 
It is possible to describe similar asymptotically uniformly most powerful tests 
even in some nonparametric situations. Using the minimax approach we can 
consider the least favorable model in the deriving of the upper bound on 
the powers of all tests, but, of course, for special classes of intensities. This 
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approach sometimes is called semiparametric and the rate of convergence of 
alternatives is a/T. 

As before, we suppose that under hypothesis the observed point pro- 
cess = {X(, < t < T} is standard Poisson with known intensity function 
S {t) = S^> ^ and under alternative M\ it is self-exciting point process with 
intensity function 



S{t,uj) = g{t- s) dXs, 0<t<T 

J — oo 

where g (■) is now unknown function. We suppose as well that 

POO 

/ 9 it) dt < 1, (9) 

hence the process is stationary. To describe the class of local nonpara- 
metric alternatives we rewrite this intensity function as 

1 /■* 

S {t,u) = + ^ u{t-s) dXs, 0<t<T 



where the function m (■) is from the set defined below. Let us denote 
by C'l, the set of nonnegative functions bounded by the same constant and 
introduce the set 

'^r = {■) E : J u if) dt = r, supp u (■) is bounded 

Note, that for any r > and T > r"^ the condition ([9]) is fulfilled. 
Therefore, we consider the following hypotheses testing problem 

JTo : u{-) = 0, 

^1 : u{-)e'^r, r>0. 

The power function of a test (pj, depends on the function m(-) and we 
write it as 

where m = m (■) G with some r > 0. We want to apply an approach 
similar to the minimax one in the estimation theory. More precisely, we seek 
to maximize the minimal power of test on the class However, for any 
test (j)rp e we have 

inf Pt (u, (brp) < e 
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since for any T > we can take a function from equal on [0, T]. Hence 
we introduce the set 

%,N = {«(■) e % ■■ suppm(-) c [o,n]}, 

denote 

Bt (r, N, (t)j) = inf (5t {u, (pj) 
and give the following 

Definition 2. A test 0^ (■) is called locally asymptotically uniformly most 
powerful in the class ^ if for any other test (•) G a-nd any X > we 
have 

hm hm inf \Bt (r, N, 01) - Bt (r, N, 0^)1 > 0. 
Let us introduce the decision function 

Theorem 2 The test (f)j, is locally asymptotically uniformly most powerful 
in the class and for any u{-) & its power function 



(3^(u,4>^^ -^/3H = p{c>^.-r v^}, (10) 
where C ~ ^/'(0, 1). 



Proof. Let us fix a simple alternative u (■) G '^r, then the likelihood ratio 
Lrp (^^^^^^ = {u (■)) admits (under hypothesis Jifo) the representation 
(see the proof of the theorem 1) 

Zt{u{-)) = exp|^^ln(^l + -^^* u{t-s)dX,^ dXt- 



exp I {u, X^)-^l (u) + r^ {u, X^) | 



where 



A^{u,X^) = -^ [ [ u{t-s)dX,[dXt-S,dt], 

O^V -L Jo Jo 

/CO / POO \ 2 poo 

u{tf dt + sjj u{t) dtj = J u{t)'^ dt + S^ 



r^ 
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and 

£o { {u, X^)}^^ (0, l{u)), {u, X^) ^ 0. 

Moreover, these last two convergences are uniform on m (■) G Ur,N, < r < K 
for any K > 0. Hence the hkehhood ratio test 

V J '^{ZT{u{-))>de} 

with = exp l^;^-^! (u) — I (m) /2| is the most powerful in the class for 
any two simple hypotheses and its power function 

(3 {u, 0^) ^ P |c > - I (uf } , C ~ (0, 1) . 

It is easy to see that 

inf I (u) = 5, + ^ 



because 



r^ = { I u{t) dt] <N / u {tf dt 



with equality on the least favorable alternative u* (t) — (r/N) X^q<«^-^- 
Hence 

inf P |C > - I W'^'l = P IC > - ry/S, + N'^X . 
Now we study the power function of the test 0^,. Let us denote 

Ut^ [ u{t-s)dXs, TTt^Xt-S.t, 

Jo 

then 

At {u, X^) = ^[ut dvr,, 5^ {X^) = dvr, 

and 

EoAt {u, X^) = 0, EoAt {u, X^Y = \ {u) , Eo6^ {X^) = 0, 
EoSr^ (X^) ' = 1 Eo {6 J, {X^) At {u, X^)) ^ r ^/^. 

Hence, under hypothesis J^, we have 

Co {At {u, X^) , 6^ {X^) } ^ A/" (0, R) 
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where R is covariance matrix of the vector (A^, 6t) described above. There- 
fore (t)j, E J^s: and using Le Cam's Third Lemma (van der Vaart (1998)) we 
obtain that under alternative u{-) e 

6^ (X^) =^ M (r^,, l) . 

For the power function we have 

(5 (u, 0y) ^ P |C > - r V:^} . 

It can be shown that this convergence is uniform over m (■) G '^r,N-, < r < K 
for any K > and this proves the theorem. 



6 Simulations 

The main results (Theorems 1 and 2) of this work are asymptotic in nature 
and it is interesting to see the properties of the tests for the moderate values 
of T. This can be done, say, by Monte-Carlo simulations. 

6.1 Parametric alternative 

To illustrate Theorem 1 we take = 1 and h (t) = |e^*/^ (see Example 1). 
This yields 

>5(ii,i,a;) = 1 + ^ ^e-^*-**)/^ u>0, 0<t<T. 



2^ u<t 



In this case 



0<tj<Tti<tj vol y 0<tj<T 

where ti are the events of the observed process, and the test 0^ given by 

^T-^T (^^)=X{^^(x-)>.,} 

is locally asymptotically uniformly most powerful in the class J^. 

In Figure 1 we represent the size of the test 0^°^ as a function of T e 
[0, 1000]. This size is given by 

a{T)=Po{AT{X^)>zom}, 1 < T < 1000 
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and is obtained by simulating M = 10'' trajectories on [0, T] of Poisson 
process of constant intensity S {t, co) — 1 and calculating empirical frequency 
of accepting the alternative hypothesis. 

In Figure 2 we represent the power function of the test 0^°^ given by 

(3^ (u, = P„ {At {X^) > Z0.05} , 0<u<5 

for T — 100, 300 and 1000, as well as the limiting (Gaussian) power function 
given by 

, ~i 1 

P{u)^f\C> zom-uV5/2\^^= e'"'/^ dv, < u < 5. 

The function is obtained by simulating (for each value of u) M = 10^ 
trajectories on [0, T] of self-exciting process of intensity S {u, t. u) and calcu- 
lating empirical frequency of accepting the alternative hypothesis. 




100 200 300 400 500 600 700 800 900 1000 1 2 3 4 



Fig. 1: Test size Fig. 2: Test power 

Now let us consider the 0^ given by 

where the threshold z is chosen so that this test is of exact size s. The choice 
of this threshold z as a function of £ e [0, 0.25] is shown in Figures 3 and 4 for 

T = 100, 300 and 1000, as well as the Gaussian threshold z^. The values of z 
are obtained by simulating M = 10^ trajectories on [0,T] of Poisson process 
of constant intensity S {t,u!) = 1 and calculating empirical 1 — s quantiles of 
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0.5 - 



0.05 0.1 0.15 0.2 0.25 

Fig. 3: Threshold choice 






0.14 0.15 



Fig. 4: Threshold choice (zoom) 



For example to obtain test of exact size 0.05 one needs take z ~ 1.78 for 
T = 100 {z ~ 1.74 for T = 300, z ~ 1.70 for T = 1000) against z^ ~ 1.64 for 
Gaussian case. 



6.2 Nonparametric alternative 

To illustrate the nonparametric alternatives we take intensity functions cor- 
responding to 5** = 1 and u (t) = (r/N) X{o<t<Ar}' 



where ti are the events of the observed process. This choice of (•) allows us 
to compare the power function of our locally asymptotically uniformly most 
powerful test 

with the asymptotic power 



oo 



Ze~r 



of Neyman- Pearson test for the least favorable alternatives. 

Note that under J^, Xt is Poisson random variable with parameter T, 

therefore the size of the test (^|^, as well as the threshold giving a test of exact 
size £, can be calculated directly (without resort to Monte-Carlo simulations). 
We represent the power function of the test 0^°^ given by 

/3^(m>t) =P«{Xt>zo.o5Vt + t}, 0<r<5 
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for T = 100, 300 and 1000 as well as the limiting (Gaussian) function /3 (r), 
< r < 5. In Figures 5 and 6 we take N — 5 and in N — 50 respec- 
tively. The function is obtained by simulating (for each value of r and 
A^) M = 10^ trajectories on [O.T] of self-exciting process of intensity S {t,uj) 
and calculating empirical frequency of accepting the alternative hypothesis. 

We see that if 1 -C iV -C T, then the power function converge to the 
hmiting function (for example, if = 50 and T — 1000, the power function 
almost coincides with the hmiting one). If N and T are of the same order 
(for example, if = 50 and T = 100) then the power function of the test 
can be essentially smaller. This example confirms the importance of use of 
functions with bounded support and of the order of limits in Definition 2. 




Fig. 5: Test power (A^ = 5) Fig. 6: Test power {N = 50) 



7 Discussions 

The constructed tests are asymptotically optimal for parametric (Section 3) 
and nonparametric (Section 5) alternatives. It seems that these are just the 
first results in this field and it is interesting to develope the construction of the 
asymptotically optimal tests for wider classes of alternatives. Particularly, it 
is intersting to study smooth alternatives like 

/•oo 

J^i : / u^''^ {tf dt > r, 
Jo 

where r > 0. Note that the test 0r is no more uniformly consistent in this 
situation. 
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